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ABSTRACT 


The results for an experimental study of a one-wavelength MHD induction 
generator operating on a liquid flow are presented. First the design philo- 
sophy and the experimental generator design are summarized, including a 
description of the flow loop and instrumentation. Next a Fourier series 
method of treating the fact that the magnetic flux density produced by the 
stator is not a pure traveling sinusoid is described and some results 
summarized. This approach appears to be of interest after revisions are 
made, but the initial results are not accurate. Finally, some of the ex- 
perimental data is summarized for various methods of excitation. 
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SECTION I 


INTRODUCTION 


The successful development of the MHD induction generator depends on 
the selection of a compromise between electrical and viscous losses such 
that an acceptable efficiency is obtained. Thus the basic objective of 
the program is to gain further understanding of the generator and of the 
associated losses. Particular emphasis has been placed on the divergence 
between the anticipated and measured performances at low slips. The philo- 
sophy was to design a generator and test facility with the maximum amount 
of flexibility and instrumentation. Steady operation was desired, and 
obtained for periods of up to a day. 

A one-wavelength-long generator geometry was selected as this appears 
to be the only feasible choice except for the largest sizes. This natural- 
ly lead to the primary emphasis being on the study of end losses and com- 
pensating methods. Provision has been made for removable compensating 
poles, grading of the traveling-field winding, and the use of both vaned 
and unvaned channels Outside of the traveling-field region. This allows 
the comparison of all cases while utilizing the same traveling-field 
structure. 

Over 200 data points have been measured for the experimental generator 
under a variety of operating conditions with the unvaned channel. These 
ciiicl ude: 

1. Uncompensated, connected to the power line. 

2. Uncompensated, connected to the power line, end coils shorted to 
reduce the effect of the open magnetic circuit. 

3. Uncompensated and self-excited for both cases above over a fre- 
quency range of about 35 to 60 hertz. 

A. Uncompensated, two-wavelength-long field connection. 

5. Compensating poles excited, connected to the power line. 

6. Compensating poles excited, self-excited over a frequency range 
of 37 to 53 hertz. 

7. Harmonic measurements and tests of system accidents. 

The data runs with the compensating poles have been limited (AO points) 
due to difficulties in properly setting the compensat ing-pol e currents, 
as different criteria (maximum power, phase balance) lead to different 
operating points. 

In Section II the design philosophy, the experimental generator, and 
the facility are briefly presented. The design is described in more detail 
in a previous report [1]*, which also includes a summary of the theory. 

The initial Fourier-series model for the MHD induction generator is pre- 
sented in Section III. The experimental results are summarized in Section 
IV. 


Numbers in brackets references at the end of the report. 



SECTION II 


EXPERIMENTAL GENERATOR AND FACILITY 


2.1 DESIGN PHILOSOPHY 

The basic objective of this experimental program is to gain further 
understanding of the operation of the MHD induction generator, and in 
particular to examine the discrepancies between the expected and the 
actual performances at low slips. To this end it is desired: 

1. To build a machine whose performance can be accurately predicted. 

2. To do a careful, controlled experiment capable of steady 
operation. 

3. To design a generator capable of reaching a high enough slip, 
say s =» -0.8 as shown by previous experiments [2], so that the 
losses which cause the discrepancies at small slips are not 
significant. 

A. To build a one-wavelength generator as this appears to be best 
for space applications. The longer wavelength of the one- 
wavelength machine yields a stronger interaction between the 
fluid and the traveling field and a better power factor. More 
than one slot per phase per pole should be used to obtain a 
better traveling magnetic field and to allow operation as both 
a single - and a mul t i -wavel ength machine by reconnecting the 
coils. 

5. To include provisions for a careful study-, of the end losses. 

The machrne must be capable of operation with and without both 
compensating poles and vanes in the end regions. 

6 . To wind the exciting coils in such a fashion that the current 
in each slot can be independently controlled, allowing grading 
of the winding for end studies and balancing the phases. This 

is accomplished by winding each coil around the back of the core 
so that each coil is in only one slot, and by including a tap 
i n each co i 1 . 

7 - To build a machine that is capable of self-excited operation as 
a generator. This is not only necessary for a practical machine 
but, lacking a variable- frequency power source, it allows 
greater experimental flexibility. It should also be capable of 
operation from a 60 -hertz three-phase power line. 

8. To include provision for operation with a two-phase (liquid and 
vapor) mixture as this is expected for practical power systems. 

These above points formed the basic framework for this generator design. 

To achieve the above major points, and to obtain the maximum amount 
of information from a single traveling-field structure or stator, a constant 
afield and fluid-velocity generator was chosen as this has more experi- 
mental flexibility than the variable-f luid-veloci ty generator with either 
constant or varying field velocity. More, important, the variable-velocity 
machine would not yield the controlled experiment desired because of the 
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interaction between the end and variable-velocity effects. It was felt 
that, the major areas where information was required were end losses 
and operation with two-phase working fluids, and that it was preferrable 
to concentrate on points 1, 5, and 8 and leave variable field velocity for 
later experiments. 

2.2- experimental GENERATOR DESIGN 

The main constraint on the size of the generator is the flow loop, 
and in particular the rating of the pump, 250 gallons per minute at 82 
psf. The generator is designed to use the full flow capacity of the 
pump only at the highest desired fluid velocity (slip), and to allow 
leeway for pressure losses in the flow loop, the pressure wasted as lost 
kinetic energy into the supply tank, and the uncertainty in the actual 
pump performance. 

The electrical constraints on the experimental generator are very 
loose. Three-phase, 60 hertz, ^80 volt power is available with more 
than adequate capacity. The generator is driven through three single- 
phase variacs for control, Capac i tors i n pa ra 1 1 el with the generator . 
are used to reduce the reactive power load on the variacs, and are 
required for self-excited operation. Operation at other frequencies is 
possible sel f -excited . 

The practical constraints for obtaining an acceptable experimental 
performance are more subtle. Some of these were mentioned in Section 2.1. 
Further constraints are: 

The viscous pressure drop must represent a small fraction of the 
total available pressure in order to model a practical generator. 
This is required both for acceptable efficiency and for a large 
field-fluid interaction where the electromagnetic force dominates 
over the viscous force. Otherwise a sharp Hartmann profile 
(or the equivalent for turbulent flow) will not be attained,' and 
the ohmic losses in the boundary layers will be unreasonably 
large [ 3 ]. 

The generator height parallel to the magnetic field, 2a^, should 
be small t6 reduce the exciting ampere turns and power required, 
but large to reduce the viscous loss. Also, has to be larger 
than the wall and/or insulation thickness so that most of the 
field volume performs useful work, inproving the power density 
and power factor. Finally, 2Tra^/A <•? 1, si i t-channel model, for 
an acceptable power density, wherex is the wavelength. 

The wavelength should be long fora strong field-fluid interaction, 
but short to reduce the viscous loss. The wavelength determines 
the room available for the exciting winding, and thus the number 
ampere turns and slots per phase per pole. 

The generator width c should be large so that the effects of the 
cooper side walls and the current loops closing through them 
are small, but increasing the width decreases the fluid velocity 
and thus the wavelength. 
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5 . 


The stator desJgn should be such as to minimize the effect of slot 
harmonics, which appear to have been important in previous 
experiments [4]. The design should allow their effects 
to be observed experimentally. 

6. Conducting channel walls should be avoided because their 
ohmic loss will be proportionally very large in a small 
generator (smaller ratio of generator volume to surface 
area). Also, the walls will be more sensitive to slot 
harmonics, which decrease away from the edge of the core. 

In order to avoid wall loss dominating the generator 
performance, hiding the very losses of interest in this 
experiment, and interfering with self-excited operation; 
it was decided to use a non-conducting fluid channel except 
for copper side bars to reduce the edge effects. 

As can be seen from the above, the design of the experimental generator 
is not a trivial exercise because many of the constraints are contra- 
• dictory . 

2.2.1 Channel Design and Dimensions 

A fluid channel is required to protect the stator and reduce the' 
danger of leakage but, as mentioned, conducting channel walls should 
be avoided. Thus, it was decided to use a fiberglass channel because 
fiberglass is flexible and has a higher temperature rating (compared 
to lucite). The minimum wall thickness a], 1/16 inch to avoid the 
danger of breakage, is chosen to minimize the required exciting 
ampere turns. It is assumed that the iron core is perfectly smooth 
and is used to mechanically back up the channel and absorb all of the 
pressure force. (The channel v/as tested at the maximum working 
pressure of 100 psi before its use with NaK) . Incidently, the core 
will probably be similarily used for structural backup in a power 
system. 

In spite of the insulating top and bottom channel walls, con- 
ducting side walls are desired to facilitate the closing of the 
fluid current loops and thus increase the power density. This is 
accomplished by using cooper side walls, as shown in Fig. 1. 

The copper bars are clamped so that they cannot be pushed out by the 
static pressure force, and have a lip which rests ajainst the core so. 
that they cannot be pulled into the channel by the J x ^ force or the 
dynamic pressure force. Leakage between the copper and the fiberglass 
is prevented by inserting the fiberglass into a slot in the copper. 

The copper extends into the region between the two cores, in spite 
of slightly increased ohmic losses, to provide additional sealing 
and the mechanical strength to support- the cores when they are clamped 
together. 

The first dimension chosen was the height a^ = 0.25 cm, compared 
to the wall thickness of 0.159 cm. The choices for X and , the 
maximum attainable slip, are independent but together fix the other 
parameters and fhe generator performance. The relevant equations 
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are v = (1 - s )fX, and c = Q /2a.Vn^j^ where f is 60 hertz and 
0 Ts^the miximum^flow rate of tRe^pump. After several trials, the 
^fues A = 0.2) m and s = -0.8 were selected, giving c = 0.1393 m. (In 
the actual generator c''^?urned out to be about 0.1405 rn because it was not 
feasible to control the width of the iron core precisely.) 

The estimated power and pressure drops for this design are shown in 
Table I, where the equations for the slit-channel induction machine ^ere 
used [3, 5 ] with the design peak magnetic flux density B = 0.4 wb/m and 
the appropriate NaK' properties of 0 ^ = 2.38 x 10 mhos/m, p = 84/ kg/m^, and 
n = 4.68 X 10“^ kg/m-sec. Only the' travel i ng-f iel d region is considered, 
perfect compensation for the ends is assumed, and the edge effects are 
ignored. In the table P is the power supplied by the fluid to the exciting 
winding (electrical output power when coil and core losses are ignored), 
is the mechanical input power, Ap is the electromagnetic pressure difference, 
and R is the Reynolds number. TBe viscous pressure drop Ap is calculated 
using^the friction factor for turbulent flow without a magnetic field [6] 
as the case with a magnetic field is not clearly documented [ 3 ]. 

The correction required for the non-flat velocity profile is shown in 

Table II. The Hartmann number M based on the rms magnetic flux density is 

50 . 4 , which is on the borderline of being too small. To see this, the ratio 

of (S*, the distance the channel wall would have to be moved in to maintain 

the same volume flow rate if the velocity was constant at the maximum (center 

line)value, to a, is calculated using an approximate turbulent boundary- 

layer solution [ 3 ]. Since for a Hartmann profile 6*/a^ = 1/M, this is used 

tocalculate an effective M, M^^^, as far as the fluid ohmic losses due to 

the non-flat profile are concerned. The effect of this is to multiply P 

by the profile factor F , which takes the extra ohmic losses into accoun?, 

but to leave P unchanged [ 3 ]. The increase in P is 7% at s , with larger 
• s , m ni 0 X 

increases as s approaches zero. 

One other point is to estimate the effect of the copper side walls or 
the finite generator width on the performance. The peak current in the 
side walls, assuming that all of the fluid current flows in them and that 
their presence does not change the fluid current from the infinite-width 
value, is 2210 amps for s = -0.8, giving an ohmic loss for both sides of 465 
watts which is limited by the skin depth of the copper. The ratio of this 
ohmic loss to that in the fluid, neglecting the profile factor, is O.I 63 
for this generator independent of the slip, so that the power density in the 
generator will be decreased. This is an upper limit, in practice the fluid 
current constraints should reduce substantially the ohmic loss in the copper. 

2.2.2, Traveling-Field Structure 

The purpose of the traveling-field structure or stator is to produce 
the required magnetic field with the minimum of distortion from a pure 
traveling sinusoid. The effect of distortion is to generate spatial 
harmonics with a shorter wavelength and slower field velocity, and conse- 
quently to reduce the efficiency due to the increased ohmic losses in the 
fluid. Distort ion- is created by many sources, the most important being 
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TABLE 1 . Powers 
Parameter 

and pressure 
Un i ts 

drops for the 

travel i ng-wave 

reg ion. 

s 


-0.8 

-0.5 

-0.3 

V 

m/sec 

22.7 

18.9 

16.4 

B 

wb/m^ 

0.4 

0.4 

0.4 

P 

s 

watts 

3570 

2230 

13A0 

P 

m 

watts 

6420 

3350 

1740 

Ap 

^em 

psi 

59.1 

36.9 

22.2 

R 

e 


3.96x10^ 

3.30x10^ 

2.86x10^ 

Ap 

f'v 

Psi 

9.4 

6.8 

5.2 

^'^total 

psi 

68.5 

A3. 7 

27.4 


TABLE 1 1 . 
Parameter 

Veloc i ty prof i 1 e 

correct ions , M = 

50.4. 

s 

-0.8 

-0.5 

-0.3 

R 

e 

3.96x10^ 

3.30x10^ 

2.86x10^ 

6*/a^ 

0.062 

0.056 

0.05 


16.1 

17.9 

20 

F 

m 

1.069 

1.084 

1 .108 
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the slots in the iron structure for the v;inding, the use of a finite and 
small number of phases with pulsed rather than sinusoidal turns densities 
and saturation of the iron. To minimize these effects, the following 
design ideas were adopted: 

1. The two stator cores were designed so that they can be operated 
with tooth opposite tooth or with tooth opposite slot in order 
to reduce slot harmonic effects and study them experimentally. 
This is accomplished by making the two end teeth of different 
widths so that reversing one core relative to the other switches 
cases. The core lamination dimensions are shown in Fig. 2, and 
the experimental fields produced by the two cases are compared 
in Reference 1. Also, short straight sections are provided at 
the end of each slot to allow the possible later insertion of 
highly-permeabi e slot plugs which would smooth out the air-gap 
magnetic field at the price of increased leakage flux. 

2. More than one slot per phase per pole or else more than three 
phases, either of which requires more slots per wavelength, 
were used, jnd this results in a smoother field as explained 
previous ly [ I ] . 

3. The width of the slot (or tooth) was tapered so that more iron 
area is avai 1 able where needed to handle the leakage flux with- 
out saturation. 

The traveling-field-structure design is only summarized here, as it 
is presented in detail in Reference I. In the experimental generator 
there was only room for 12 slots. The number of ampere turns was cal- 
culated to be 1230 per slot at s = -0.8, the maximum case, including 
the Carter coefficient and a small allowance for iron reluctance. Each 
coil consists of AS turns made out of three §\h wires in parallel in the 
slot butt-soldered to a rectangular bar of about twice the copper area 
outside the slot. The coil loss was estimated to be kkO watts for this 
case using an assumed effective copper area of O .7 and a larger copper 
area outside of the slot. (This is low because the effective area was 
off, and will be increased to about 700 watts.) The core loss for the 
Armco electrical steel M-lA was estimated to be 32 watts including the 
lamination packing factor. 

One important point in the stator design is to build in adequate 
monitoring devices. Two types are required - thermocouples to measure 
the coil temperatures at various locations in the slots and in the loop 
around the back of the core, and coils to monitor the field. There are 
19 two-turn magnetic-field monitoring coils on each core, located as 
shown in Fig. 2. The locations were chosen to observe the field in the 
core above the slots, which includes the compensating flux (see the next 
section); the flux through some of the teeth, including the leakage flux; 
the flux through some of the slots; and the flux entering the air gap. 
They will be useful in evaluating the stator design and in monitoring 
the operation v/ith fluid present in the channel since then the field 
cannot be directly measured. 
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2.2.3 Compensating Poles 


In the regions where the fluid enters and leaves the traveling magnetic 
field undesired voltages and currents are induced in the fluid, and these 
cause increased ohmic losses. In order to set up the proper boundary 
conditions so that the traveling-field region acts like a segment out of 
an i nf i ni tely-long machine, Elliot proposed the use of compensation poles 

[7] - The peak magnetic flux density B and length L must satisfy 

B L = BX/2tt assumi ng a slit-channel machine. Arb i trar i ly choos i ng L 
tS Se one inch gives B =0.525 wb/m^. The channel wall thickness urider 
the compensating polesSvas increased to l/8-inch to give a stronger 
structure because the poles are removable. Then, 2380 ampere turns is 
required for each of the, four pole structures ,. and 85 turns with a tap 
at 70 turns are on each pole. Excitation is by means of the phase voltages 
rather than the phase currents because this tends to be self-compensating 
for the phase shift required as the fluid reaction field varies. (This 
is discussed in Reference I.) 

Insulating vanes in the channel under the compensating poles are 
required to reduce -the ohmic loss, but at the penalty of increased 
viscous loss. It is necessary to calculate (of estimate) both losses as 
a function of the number of vanes and optimize. The viscous loss is 
easily calculated using the available hydrodynamic equations, but the 
ohmic loss can only be crudely estimated. Using the results of Moszynski 

[8] to estimate the ohmic loss, a compromise of six vanes was selected. 

l.l.k Discuss ion 

The results of the previous calculations are summarized in Table III 
to give a tentative idea of the performance to be expected at three 
values of slip. All of the losses are considered, except that the field 
harmonics are not included and perfect end compensation is assumed. 

The effect of the copper side bars is to reduce the fluid current and 
thus the power density. This effect is estimated by assuming that the 
induced electric field or voltage in the fluid stays constant since the 
magnetic flux density is assumed to be constant. (This means that the 
sum of the ohmic losses in the fluid plus the side bars decreases as the 
bar conductivity decreases.) The ratio of the ohmic loss in the side 
bars to that in the fluid neglecting the velocity profile correction is 
defined as B, where B = O.I 63 for this generator from Section 2.2.1. 

Then the ratio of the equivalent resistances of the two regions is also 
B, and the fluid current is multiplied by 1/(I+B) for the same induced 
^o^tage. Finally, the generated power in the fluid, which depends on 
E*J, is P /(I+B); where P is the power without the side bar correction, 
Table I. fhe mechanical input power modified for the side bar cor- 
rect! on is P ' . 

The separate power losses are also included in Table III. The ohmic 
loss in the exciting winding is based on the measured DC coil resistances. 
The reduction in the fluid currents due to the bars also appears here 
because the reaction field and thus the required ampere turns are reduced. 
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Table 111 . Summary of theoretical generator performance. 


Paramenter 

Units 




s 


o 

CO 

-0.5 

-0.3 

Ii/(i+e)]p^ 

watts 

3070 

1920 

1 150 

P 

e 

watts 

6 k 0 

590 

560 

p 

watts 

90 

90 " • 

90 

core 

P 

CP 

watts 

230 

230 

230 

P ^ 
out 

watts 

2110 

1010 

270 

p'/p 


1.67 

1 ,42 

1 .26 

ms 

P* 

m 

watts 

5980 

3160 

1680 

Ap 

“em 

.PS' 

55.0 

34.8 

21.4 

APv 

Psi 

9 .i» 

6.8 

5.2 

^^cp 

ps i - 

10.8 

7,5 

5.6 

^Ptot 

psi 

75.2 

49.1 

32.2 

'“in 

watts 

8190 

4460 

2530 

e 


0.26 

0.23 

0.11 

g 
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This is taken into account by multiplying by l/(i + g) , and amounts 

to decreases in excitation loss of 8^, and 3^ for s = -0.8, -0.5» 

and -0.3. The core loss, P , is essentially independent of slip 
or the side or the side bar^effect except for the slight change in 
the leakage flux which is ignored. The excitation power for the 
compensating poles, P , depends only on the flux density which is 
independent of the sl*Tp. The net electrical output power is 


Similarly, the mechanical pressure drops are taken from the 
previous results. The electromagnetic term Ap comes from P', and Ap 


is given in Table I. The pressure drop for the'^compensat i ng poles, 
Ap , is estimated by assuming that the electrical drop is equal to 
the*^viscous drop since the calculation for the electrical pressure 
difference was very rough. The net result is estimated efficiencies 

of 0.26 at s = -0.8, 0 . 23 , at -0-5. and 0.11 at -O. 3 . 


2.3 THE FLOW LOOP AND INSTRUMENTATION 


The existing NaK-N^ flow loop at the Argonne National Laboratory 
with the addition of a new extension intended specifically for the 
induction generator is used for these experiments. The pump, heat 
exchanger, supply tank, separation system, and part of the control 
system are utilized from the existing loop. The flow path is from the 
supply tank to the pump, the heat exchanger, then to the extension 
and generator, and back to the storage tank. For the two-phase ex- 
periments the supply tank also serves as the separator. There is 
an air-operated shunt or bypass valve between the pump outlet and the 
supply tank, but this valve is always closed for tests. 

The pump, rated at 250 gpm at 220 ft. head, has a combined 
rotor-impeller assembly with a completely sealed or canned stator 
winding. A small portion of the pumped NaK is allowed to recirculate 
through the rotor cavity to cool the motor and lubricate the 
aluminum oxide hearings. This construction leads to a compact pump 
with no possiblity of leakage, but since it dumps nearly all of the 
motor ohmic loss into the NaK it may also cause an excessive temperature 
rise in the NaK. To maintain the NaK temperature below about 125°F, 
required f6r the fiberglass tests section, a stainless-steel 
finned-tube NaK-to-air heat exchanger is' used. The cooling capacity 
is approximately 12 kilowatts at an outside air temperature of 75°F, 
a NaK inlet temperature of 150°F, and a flow rate of 100 gpm; but 
this is less than the power dumped into the NaK by the pump. 

The extension for the i nduct ion" generator tests contains a 
venturi coupled with a differential pressure transducer for measuring 
the NaK flow rate, a three- inch air-operated valve to control the 
pressure across the test section, and three- inch ball valves to y 

isolate the extension from the rest of the loop and from the test 
region so that system changes can be easily made. The test region includes 
the mixer for adding nitrogen to the NaK flow and the fiberglass test 
section. The nitrogen supply will be obtained from a ^2,500 ft^ 



capacity trailer. The nitrogen is fed into the mixer through four check 
valves to prevent flow of the NaK back into the nitrogen system. The 
mixer consists of an 11-inch length of three-inch Schedule lOS stainless 
steel pipe with 900 1/8-inch diamerer holes enclosed by a concentric 
length of five-inch Schedule kO pipe into which the nitrogen is 
introduced. The nitrogen flow will be measured with an orifice and a 
differential pressure transducer. It will also be necessary to measure 
the void fraction by means of a gamma-ray source and a scintillation 
photo-multiplier tube. .y 

Pressure measurements in the test section are made by means of 
gauge-diaphragm seal assemblies. A 316 stainless steel diaphragm 
is welded to the upper chamber in a two-chamber, 316 stainless steel 
pot assembly. The upper chamber is connected to a gauge and/or a 
transducer by tubing and the upper assembly is filled with instrument 
oil. The gauge accuracy is \ H% of full scale. The lower chamber 
of the diaphragm unit is connected to the test section and to a 
bleed system for purging. Five pressure taps are provided on 
one side of the channel with a spare set on the other side. The tap 
locations are at" the entrance, exit, and center of the traveling-field 
region and the outer ends of the compensating poles. 



SECTION III 


AN INITIAL FOURIER-SERIES MODEL* 


3.1 INTRODUCTION 

Several different approaches have been utilized to develop analytical 
models for the MHD induction generator. Initially most studies were based 
either on field models using a single traveling magnetic field or on the 
techniques developed previously for rotating induction motors, such as 
circle diagrams. These were limited in their treatment of slot harmonics, 
fringing fields, fluid velocity profiles C3D> etc. Slot harmonics, for 
example, can be minimized for rotating induction machines because the 
slots in the rotor can be skewed, but harmonic effects may be significant 
for linear induction machines. Recently two new approaches to linear 
induction machines have been proposed, and these can treat arbitrary cur- 
rent-excitations and stator geometries. They are: (1) Elliott USD 

analyzes the induction generator by breaking the stator winding and the 
fluid secondary into discrete pieces and then solving for the total current 
distribution and magnetic field using superposition and matrix techniques. 
(2) The Four ier-series method described below. 

3.2 EMPTY-CHANNEL MAGNETIC FLUX DENSITY 

The magnetic flux density with no fluid present is calculated in three 
steps for the tooth of one core opposite the slot of the other core. This 
is the only case studied experimentally with the fluid present. First the 
field due to the current in the coil sides in a pair of slots is determined 
and broken down into Fourier series components. A slit-channel model is 
used, that is the magnetic field is assumed to be completely transverse to 
the fluid velocity and to not vary over the transverse direction. The 
neglect of the transverse variation means that the alternate location of 
the coils and slots on either side of the air gap has no effect on the 
field. In addition, fringing at the ends of the stator (magnetic core) 
is neglected, the effect of slots in the stator iron is not considered at 
this point, and the iron is assumed to be infinitely permeable. Second, 
the total field is found by adding together the contributions from the 
currents in all 2A slots or coil sides, first for the eight per phase and 
then for the three phiases. Finally the effect of the slots in the iron is 
included and the total magnetic field appropriately modified. 

3.2.1. Magnetic Field Due to a Pair of Coil Sides 


The magnetic flux density produced by. a single coil located 
slot is found from Ampere's law 



C S 


Based in part on a. special project performed by M. A. Brown. 


in a single 


(3.1) 
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conservation of flux, 

f 

0 (^*^)da = 0, (3*2) 

4 

5 ^ . 

and the constituent relations for .the iron and the air gap, where B is the 
magnetic flux density, J the current density, and p the permeability. For 
a slit-channel model and infinitely permeable iron only the transverse 
flux density B has to be considered. With no magnetic fringing Eq. 3*2 
s impl i f ies to ^ 

j Bdx =0, (3.3) 

0 

where L is the length of the stator. 

With no fluid present the slit-channel single-coil field has no spatial 
dependence except at the coil location where the field discontinuity is 

AB=p^NI/g (3.^) 

for Nl ampere-turns per coil (N turns carrying I amperes) and an air gap 
of g. The time dependence is contained in I, and the coil is assumed to 
have negligible width. _ 

The field due to the combination of two coil sides located in slots 
a half wavelength apart is found from the jump and conservation conditions, 
Eqs. 3.^ and 3.3* The form of the field is a square wave with different 
amplitudes for the positive and negative excursions to preserve the zero- 
average value. The resulting Fourier series for the two coil sides is 

p N I “ 

B = -^—l [c^cos(5mx) + D^sin(Cmx)]; 

^ m=l 

where 

^m ~ m [>^5 ~ s in [m? (d+(n+l 2) 6)] ■, 

Dm = ^ "Icos [m^ (d+(n+l 2) 6)]] - cos [m^ (d+n6)] •, 

C = 2ti/L 


(3.5) 

(3.6) 

(3.7) 

(3.8) 


is the wave number for the series, and n is the number of the first coil side 
under consideration (0 <. n ^ 1 1 ) . Here L, the stator length, is 8.947 
inch and is made up of 24 slots spaced 6 =.0.3445 inch apart and two over- 
sized ends of length d = 0.512 inch. Note that ^ is not the same as the 
wave number k = 2i\/X for the excitation winding due to the stator con- 
struction used. 

An alternative approach to that used above is to determine the field 
due to each coil side independently. This allows the results to be easily 


14 



extended to include the compensating coil currents and methods of grading 
the excitation. 

3.2.2 Magnetic Field Due to All Coil Sides 


The exciting currents for the three-phase system are 


*A ~ IgCos(tot), ( 3 - 9 a) 


and 


I g = IjjCOS (cot-(j)) , 
I . = I cos (ojt-e) 

V/ C 


( 3 . 9 b) 

( 3 . 9 c) 


where I , 1 ^^ and I are the amplitudes of the three currents and (fi and 6 
are the^ phase angles for Ig and 1 . relative to 1 ^. For a balanced three- 
phase excitation the three amplitudes are equal, (j) = 120*^, and 0 = 2 k 0 ° . 

The best experimental approximation to a pure traveling sinusoidal 
magnetic field and the best phase relations for the exciting currents were 
obtained with a symmetric connection of the coils Cl 0 ,ll 3 . In this case 
the coil at each end and the two central coils from each stator were con- 
nected to the same phase, with the coils from the two stator cores in 
parallel. The resulting coil currents for the generator, since the upper 
and lower core slots alternate, are thus -I. in coil sides (slots) 0 and 1 
(numbering from 0 to 23 to correspond to the definition of n) , +lg in 
coils 2 through 5, ~l- in coils 6 through 9, +1^ in coils 10 and 11, and 

then the negative of the preceeding currents in^l2 through 23. 


The total magnetic flux density for all three phases without any fluid 
present or slot effects included is 


u N 

R = ° 

E 2irg 



m= 


AjQcos (m^x+ut) + cos (m^x-wt)] 

+ A2Csin(m^x+o3t) + sin (m^x-oit)] 

Ip “I 

Bj[[cos (m^x+ajt-iJ)) + cos (m^x-wt+ij)) J 
B2CS in (mCx+(ot-(j)) + s in (m 5 x~d)t+({))li| 


m L 


CXcos (m^x+ut-e) + cos (m^x-ut+e) ] 


+ C2Csin(mCx+a)t-0) + s i n (m 5 x-o)t+ 0 ) U 


>]}= 


(3.10) 


where 




(3.11a) 
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(3.11b) 


Az = -F^(0) - F^(l) + F^(IO) + F^(ll), 


5 


^1 = 

- 1 fsl"). 

n=2 

(3.11c) 

^2 = 

Z %(n). 

n=2 

(3. lid) 

S = 

* 

!%("). 

n=6 

. (3. lie) 


- Z 

n-6 ^ 

(3. Ilf) 

Fg(n) 

= sinCm5(d+(n+12)6)3 - sln[mC(d+n6)3 , 

(3.12a) 

F^(n) 

= cosCmC(d+(n+12)6)3 - cosCrn^ (d+n6)3 . . 

(3.12b) 


The final step in calculating the f iel d wi thout the fluid and slot 
effects is to include the finite width of the coil. in the above the coil 
width has been considered to be an impulse. To include the finite coil 
width a five-point average of the field given by Eq. 3.10 is taken in the 
numerical evaluation to convert the steps in the field at each slot to 
jumps with a non-infinite slope. 

3 

3.2. -i Finite-Width Slots 

The effect of the slots in the iron cores for the coil sides ..is to 
decrease the magnetic flux density under each slot. This is taken into 
account in the analysis by multiplying the field given by Eq. 3-10 by a 
modulating function due to the increased air-gap length. The slot effect 
will depend on where between the two cores the field is measured or calcu- 
lated. Here attention is confined to the center line of the channel where 
the modulating functions for the slots in both cores are equal. 

The trial modulating functions were chosen because they had approxi- 
mately the correct shape as determined by Carter ■[] 123 • The first function 
tested was a gaussian density, but a computer solution proved to be im- 
practical. Thus a sinusoidal function was used with the arguement adjusted 
so that the modulation due to each slot decreases to zero at the edge of 
the slot. Because the slots are wider than the teeth and thus overlap 
slightly, the total modulating function is always less than unity except 
outside of the slots at the ends of the stator. 
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The modulating function chosen is 

23 

f(x) = n Cl - af,(x,n)Il (3*13a) 

n=l ' 


where 


(x,n) 


cosC''^ (x-d-n6)2wD, |x-d-n6| < w 
0 , |x-d-n6 I > w 


(3.13b) 


w is the slot width, and a is the modulation amplitude. The resulting 
expression for the total empty-channel magnetic flux density is 


= f(x)B^ 


( 3 . 1 ^) 


where B^ is given by Eq. 3.10. 

The modulation amplitude was obtained experimentally by matching the 
average values of the_experimental field and the theoretical field as given 
by Eq. 3.1^ for the same exciting current amplitudes and phase angles. 

The value of oi obtained, 0.l6, corresponds to a Carter coefficient of 1.22. 
This is very close to the value of 1.20 used in the generator design for • 
the actual ratio of slot width to air-gap length CO. 

The peak magnetic flux density measured with a 20-turn i-inch-diameter 
pick-up coil is shown in Fig. 3 for ^0 amps per phase, current phase angles 
of - 127 . 3 °, and + 129 . 8 °, and no channel or fluid. Measurements were made 
every 1/16 inch. The calculated field using the measured exciting current 
amplitudes and phase angles, also shown in Fig. 3, agrees well with the 
measured field. The calculated field is shown as it would be measured by 
the same pick-up coil. The calculated field is larger hear the ends of the 
stator because of the neglected fringing, and this also has an effect away 
from the ends because of the zero-average field condition used for each 
coil. 


3.3 FLUID-FLUID INTERACTION 

The empty-channel magnetic flux density will be attenuated and shifted 
in phase by the addition of the fluid. One advantage of the Fourier series 
approach is that the effect on each component of the Fourier series |s 
already known ElU. The amplitude is multiplied by Cl + (sRj^2a^/g) ? 
and the phase is shifted by arctan (sRj^2a r/g) . The factor 2a^/g appears 

in both expressions because that is the traction of the total air gap 
occupied by the fluid. Here s = (v - v)/v is the fluid slip, v and v 
are the magnetic field and fluid vefocities^ and R^ is the magnetic Rey- 
nolds number. Both v^ and R„ are functions of the field harmonic m, and 
are given by 


v^ = o)L/2Trm 


(3.15) 

( 3 . 16 ) 


and 


• Ru = P o,v L/2Trm 
n o f s 
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PEAK MAGNETIC FLUX DENSITY (N TESLA 


where o is the electrical conductivity of the fluid. The total stator 
length, L, rather than the fundamental wavelength, A, hs the appropriate 
length. 

The magnetic flux density with the fluid present is given by 


^T " Yng S 




m I, 


AjCcos (m 5 x-ij^+tot)+cos (m 5 x-'j^-ajt)H 


+ A 2 CS i n (m^x-i^+ojt) +s i n (m£x-i()-a)t) 3 


>] 

i.r 

+ — B,Ccos (m 5 x-iJ^+wt-(j))+cos (m^x-ij^-a)t+(}i )3 
m |_ I 


+ B^Cs i n (mCxr'Jj+ut-()>)+s i n (m^x-ij^-wt+ij)) 3 


■|cjCcos (m5x-i|'+tjjt-0)+cos (m5x-iJi-a)t+6) 3 


m 


.+ C^Cs i n (mCx-i|^+a)t-8.)+s in (mCx-i/^-cot+O) 3 

(3.17) 



Fig. 3 - Measured and calculated peak magnetic flux density amplitudes for 
symmetrical wye-connected excitation. 
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where 


(3.18) 


ip = arctan(sRj^2a^/g) 

and the coefficients are given by Eqs. 3-11 and 3-12. For the experimental , 
generator, the addition of the fluid has little effect on the magnitude 
and causes a small phase. shift. 

3A power 

The powers are calculated from the Fourier series terms by using the 
power-flow formulas Cl,3ll for each magnetic-field term in the series with 
the correct slip andrmagnetic Reynolds number. The power contribution from 
each term is orthogonal, so tha*t superposition of the powers can be used. 
This is one major advantage of the Fourier series technique. The calcu- 
lated magnetic flux density with the fluid present, the net electrical 
input power, the net mechanical output power, and the ohmic loss in the 
fluid are listed for the first three harmonics (m = 1 through 3) or six 
terms of the series in Tables IV to VI for three experimental data points. 
The electrical and mechanical powers are negative for generator operation, 
and excitation and viscous losses are not included. The negat i ve ‘harmon ics 
are fields traveling in the direction opposite to the fluid. Note that the 
fundamental (+1) makes the major contribution to the totals powers, but 
that the other harmonic contributions become relatively larger at the 
smaller slip values desired for high-efficiency generator operation. As 
the slip based on the fundamental approaches zero, the fundamental powers 
approach zero but the other terms change little. For the s =-0.2A case. 
Table IV, the harmonic contribution to the ohmic fluid loss is equal to 
that for the fundamental, but the percentage contribution to the other 
powers is smaller. The harmonic contributions change relatively little with 
the fundamental slip, but their percentage contribution is less for the 
higher-slip cases listed. 

The measured and calculated net electrical powers are compared in 
Table VII for five data points, where the estimated ohmic coil ioss has 
been removed from the measured power. The differences between the,, two are 
substantial, Indicating that refinements in the calculations are needed. 

The discrepancy is largest at the smallest slip value, as expected. Note 
that core, side-bar, and velocity-profile losses are not included in the 
calculated powers. 

3.5 DISCUSSION 

Improvements in the existing Fourier series approach have not been 
fruitful because some 6f the basic assumptions are not good. Rather than 
make minor improvements it was decided to revise the procedure. The 
magnetic field due to the current in a single slot is calculated first, 
including the effects of the slots and the fringing due to the finite 
stator length on the magnetic field before applying the zero-average 
magnetic flux condition. If this is not done in the proper order the 
amplitudes obtained will be incorrect, and this is particularly important 
for slots located near the ends of the stator. A single slot is chosen 
rather than a pair because then various methods of grading the exciting 
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current may be evaluated theoretically. Superposition is used to obtain 
the total magnetic field without the fluid. The effect of the fluid on 
the magnetic field and the resulting powers are then obtained as above. 
This technique is presently under development for the experimental 
generator, but it may be applied to any linear induction machine. 


Table IV. 

Fourier Series 

1 powers 

i for 20 amps per phase, s = 

-0.24, 

and 


current phase 

angles 

of -131.5° 

and 

128.4° 



Ha rmon i c 

Magnet i c 


Net Elec. 


Net Mech. 


Fluid 


Flux 


input 


Input 


Ohml c 


Dens i ty 


Power 


Power 


Loss 

+ 1 

0.277 


-6^»3 


-795 


152 

-1 

0.019 


30 


-37 


67 

+2 

0.026 


-9 


-22 


13 

-1 

0.009 


3 


-6 


9 

+3 

0.018 


-3 


-13 


9 


0.010 


2 


-7 


9 

TOTAL 



-621 . 


-919 


298 

Table V. 

Fourier Series 

powers 

for 20 amps 

per 

phase, s = - 

•0.74, 

and cur- 


rent phase angl 

es of - 

132.7° and 

129. 

7° 



Harmon i c 

Magnet i c 


Net Elec. 


Net Mech. 


Fluid 


Flux 


1 nput 


Input 


Ohmic 


Density 


Power 


Power 


Loss 

+1 ■ 

0.235 


-l^tSO 


-2530 


1070 

-1 

0.020 


38 


-66 


104 

+2 

0.02if 


-13 


-45 


32 

-2 

0.009 


3 


-11 


14 

+3 

0.018 


-5 


-27 


22 

-3 

0.010 


3 


-14 


16 

TOTAL 



-1420 


-2760 


1340 

Table VI. 

Fourier Series 

powers 

for 20 amps per phase, s = 

0.41 , 

and 


current phase 

angles 

of -130.9° 

and 

126.7° 



Harmon i c 

Magnet i c 


Net Elec. 


Net Mech. 


Fluid 


Flux 


Input 


1 nput 


Ohmic 


Dens i ty 


Power 


Power 


Loss 

+1 

0.267 


1030 


608 


426 

-1 

0.015 


13 


-8 


21 

+2 

0.027 


-1 . 


-1 


0 

-2 

0.009 


2 


-2 


4 

+3 

0.018 


-1 


-2 


1 

-3 

0.010 


1 


-2 


3 

TOTAL 



1050 


585 


464 
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Table VII. Comparison of measured and calculated net electrical input 
powers. 


Fluid 
Veloc i ty 


Slip 


Net Electrical Input Power, Watts 
Measured Calculated 


4.4 

0.68 

2020 

1430 

8.0 

0.41 

1720 

1050 

16.9 

^ 0.24 

765 

-621 

18.6 

- 0.36 

450 

. -909 

23.8 

- 0.74 

-480 

-1420 
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SECTION IV 


EXPERIMENTAL DATA 


k.] introduction 


Concurrently with developing an improved analytical model, the experi- 
mental data is being processed for consistency and correlations. As de- 
scribed in the introduction, there are roughly 200 data points for a variety 
of excitation conditions, both with and without the compensating poles. All 
data obtained to date has beenrtransferred to IBM cards, and is available to 
others developing analytical models. (With the present emphasis on analysis, 
data exchange on the existing experiments should, be encouraged.) Each data 
point consists of 22 cards containing 1^3 to 196 numbers, depending on the 
case, for the terminal voltage and current magnitudes and phase angles, the 
powers, the 38 magnetic-field monitoring-coil voltages, pressures, etc. It 
requires at least twenty minutes for one data run, but generator operation 
is reasonably constant over the time period. 

The data sheet f<3r the data point. 20-4-28 is shown .in Fig. 4, and the 
computer printout of the data in Fig. 5. The zeros on card 12 are where the 
data for the self-excited operations is entered. The case shown is for 
operation connected to the povjer line and with the compensating poles present. 

Checking the data for consistency involves utilizing the built-in re- 
dundency among some of the data, double measurements, and plotting data to 
eliminate obvious errors. Some errors have been detected, but the overall 
results have been very good. As an example, particular difficulty during 
operation was experienced with the phasemeter. This is a key instrument, 
not needed for dc, because without phase data.it is impossible to reconstruct 
the magnetic field from the exciting currents or from the magnetic-field 
monitoring-coil measurements. It has been necessary to correct some data for 
a 180*^ inconsistency of the instrument, but otherwise most results have 
checked within i2° between the readings, the built-in redundency, and the 
power calculations except for a few obvious errors or cases where the 
instrument would not yeil'd a reading. From the total set of data points 
only a few have been rejected as incomplete or inconsistent, plus some 
early cases where the measurements indicated the strong possibility of a 
two-phase mixture in the generator Cl ill- 

The measured and calculated phase angles for the data point 20-4-28 are 
shown in Table VIII. The first three angles are the voltage-current angles 
for the three phases, the others are the angles between the line voltages. 

The calculated angles use either the redundency or the power factors. 

4.2 OPERATION CONNECTED TO THE POWER LINE 

Data points were measured for various terminal conditions and excitation / 
currents between 20 and 40 amps per phase. Attention is confined here to 
the results for 20 amps because higher fluid velocities could be attained. 

Also pump operation at higher currents was limited due to excessive vibra- 
tion and/or cavitation. The powers and efficiency are increased for higher 
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Fig. A. Data Sheet for Data Point 20-4-28. 
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Computer Printout of Data Point 20-A-28. 



phase currents, but the qualitative behavior is the same. 

The net electrical input power to the fluid (gross electrical input 
power minus excitation loss) for the first 79 data points at 20 amps per 
phase without the compensating poles is plotted versus fluid velocity in 
Fig. 6. The first 18 data points (triangles) were measured before adding 
more NaK to the flow system. The evidence indicates that for these cases 
the working fluid was a mixture of liquid and vapor, so that the venturi 
used to measure the fluid velocity read low. (The static pressure is much 
less in the generator than in the venturi, so that the mixture expands in 
the generator.) The next 46 points (squares) were measured after adding 
more NaK. The last 15 data points (circles) correspond to a change in the 
technique used to measure the fluid velocity which is more accurate at low 
velocities. Each group is quite consistent within itself for an experi- 
ment. The data is now being corrected for the velocity measurements. In 
this fashion data that would otherwise be lost is utilized. 

At zero slip, 12.6 m/sec, there is a significant discrepancy from zero 
net power. This is probably due to end effects, the copper side bars, the 
fluid velocity profile, the slot harmonics, and the poor current phase angles 
(around 10° off from 120°). The improved computer program will eventually 
help to resolve the source of this discrepancy. (See Table VII.) 

. Data points were also measured with the end coils on each stator shorted 
to reduce the magnetic flux through the end teeth, or in a sense to allow 
the generator to provide its own compensation. This case does not appear to 
be attractive because the slope of the power versus velocity curve is reduced 
from that for the previous case, or there is less useable power. This was 
also true for self-excited operation. 

The generator was also operated with the compensating poles installed 
for a variety of terminal conditions. Data was taken first with the com- 
pensating windings open circuited and short circuited, as this will be 
useful in determining the characteristics of the poles and their interaction 
with the main winding. It was not possible to pick an optimum setting for 
the magnitude and phase of the compensat ing-pol e currents, as adjusting them 
for maximum power lead to inbalances between the three main current phase 
angles, while balancing the phase angles unbalanced some other quantity. 
Further analysis of the data is required to understand the operation of the 
compensating poles with this generator and the optimum setting of their 
excitation. 

Table VIII. Angle checks for data point 20-4-28. 

Measured Calculated 

From other Angles From Power 


74.2 

'■ 

73.2 

76.8 

76.8 

76.8 

73.6 

73.7 

72.5 

-117.4 

-117.4 


118.9 

119.3 


- 123.6 

-123.3 
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i».3 SELF-EXCITED OPERATION 


The generator will operate in the self-excited mode, no connection to 
the power line, with capacitance connected across the terminals to supply 
the required reactive power. The generator then operates as an oscillator 
producing net output power at the frequency determined by the generator in- 
ductance and the external capacitance. The main points noted were: 

1. Operation is very stable, and there is no trouble with fluctuations 

2. The generator starts, even after several days, without any external 
trigger. All that is required is external capacitance and a fluid 
velocity sufficiently above the field velocity determined by the 
capacitance. 

3 . The frequency is controlled by the external capacitance, with a fre- 
quency change from no load to full load at a given capacitance 
value of only a few tenths of a hertz. • 

k. The generator is short-circuit safe because it loses its excitation, 
and restarts ‘ upon removal of the short. The variation in output 
voltage with load may be decreased by adding capacitance in series 
with the loact, although this was not tested. 

5 . Operation appears to be better self-excited, probably because the 
generator adjusts the current phase angles to minimize the losses 
or optimizes the performance within the terminal constraints. The 
generator appears able to do this better than an external operator 
or control system.. 

The external capacitance on the three phases could be adjusted to bal- 
ance the phase currents, but the powers for the phases were not balanced, as 
occurred for normal line operation, and optimizing the output power required 
different load resistances for the three phases. The fluid velocity or slip 
varied with the generator load but was almost independent of the control 
valve, provided that self-excited operation could be maintained. The start 
of self-excited operation is not instantaneous, but requires a few seconds 
before anything happens. 


The known external capacitance C of the generator and the frequency f 
can be used to calculate the equivalent inductance L of the three phase in 
parallel from the conventional oscillator equation 


2it/LC 


(i».l) 


The inductance is about A. 7 millihenries without the compensating poles, and 
is constant over the 35 to 60 hertz frequency range tested. This is about 
twice the inductance calculated from the stored magnetic energy in the air 
gap between the stators. The difference must be due mainly to the leakage 
flux. 


According to Eq. A.l the plot of log C versus log f should be a straight 
line. This is shown in Fig. 7 for operation without the compensating poles. 
The slight spread over frequency at a single capacitance value is due to the 
slight change in frequency with load or output power and fluid velocity. For 
example, the spread around log f = 3.92 is from 50.4 to 50.9 hertz, or one 
percent. 
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Fig. 7. 


Log (capacitance x 10^) versus Log (frequency) for self-excited 
data points without compensating poles. 
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Seven data points have been obtained with the compensating poles in- 
stalled and excited. It was not possible to draw any conclusions from these 
because of the difficulty in defining the optimum settings, as discussed in 
Section ^.2 for operation connected to the power line. More study of the 
compensating pole operation is required. 



SECTION V 


CONCLUSIONS AND RECOMMENDATIONS 

A considerable amount of experimental data has been accumulated for 
this experimental MHD induction generator. The data is internally con- 
sistent, but interpretation of it is restricted due to the limitations of 
the existing analytical model. 

The Fourier series analysis of the generator described in Section III 
appears to be promising once the required changes in the procedure are 
included. The available method is interesting but not accurate enough for 
practical use. 

This research program on the MHD induction generator should be extended 
in two parallel paths. First the Fourier series model must be revised, and 
then compared against the existing data. An accurate analysis technique is 
needed to interpret and fully utilize the existing data. Second, there are 
several further experiments to perform with the existing stator before the 
design and construction of a new stator. Additional study of end loss and 
compensation methods is needed, particularly further data with the compen- 
sating poles excited. Also tests have never been performed with the channel 
with vanes under the compensating poles to suppress end currents due to 
leaks in the channel, and these tests are very important. Finally, a series 
of studies on a two-phase working fluid using the same stator is an impor- 
tant corollary to the above program as this will allow the comparison of the 
single- and two-phase results without the complication of a different 
generator structure. 
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SECTION VI 


NEW TECHNOLOGY 

No new technology items have been identified i n.'connect ion with this 
program. 
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